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In this paper it is demonstrated that based on two-dimensional distributions in semi-inclusive 
deep inelastic scattering (SIDIS) data, obtained recently by the HERMES experiment at DESY on dif- 
ferent nuclei, which contains data for charged pions produced in z slices as a function of v, and in 
a V slices as a function of z, it it is possible to parametrise a ratio of multiplicities on nucleus and 
deuterium (per nucleon) in a form of a function of a single variable r, which has the physical 
meaning of the formation time of hadron. We call this effect r scaling, r is a function of two variables 
V and z. It is also shown that i^Jjj can be presented in a form of a linear polynomial of r, an + ran, 
where parameters an and an do not depend on v and z. 
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I. INTRODUCTION 

Hadronization is the process through which par- 
tons, created in an elementary interaction, turn into the 
hadrons. Experimentally the hadronization process in 
free space (vacuum) has been studied extensively in 
e+e" annihilation and in semi-inclusive lepton-proton 
deep inelastic scattering (DIS). As a result, the spec- 
tra of hadrons produced and their krnematical depen- 
dences are rather well known. However, little is known 
about the space-time evolution of the process, because 
the products of this process can only be observed in a 
detector that is separated from the reaction point by a 
macroscopic distance. It is worth to mention that ac- 
cording to theoretical estimates the hadronization pro- 
cess occurs over length scales that vary from less than 
a femtometer to several tens of femtometers. The nu- 
clear medium can serve as a detector located directly at 
the place where microscopic interaction happens. Con- 
sequently, leptoproduction of hadrons on atomic nuclei 
provides a way to investigate the space-time picture of 
the hadronization process. The semi-inclusive DIS of 
leptons on nuclear targets is widely used for the study 
of this process 1 1, 2, It is most effective to observe at 
moderate energies of the virtual photon, when the for- 
mation time of the hadron is comparable with the nu- 
clear radius. Such possibility provides HERMES experi- 
ment at DESY, which uses electron (positron) beam with 
energy 21 .hGeV and fixed nuclear targets. 

The most convenient observable measured experi- 
mentally for this process is the nuclear attenuation ratio, 
which is a ratio of multiplicities on nucleus and deu- 
terium (per nucleon) for a given hadron. We shall de- 
note it as For a more profound understanding of 
the hadronization mechanism, it is important to find a 
variable which allows to present this observable in the 
most simple functional form. 

Usually it is supposed, that FC^i is a function of two 
variables v and z, which are the energy of photon and 
the fraction of this energy carried by the final hadron 



with energy (z = Eh/v), respectively^. 

In our preceding work we performed a fit for the 
evidence that the formation time r is the best variable 
for i.e. that it can be parametrized as a function of 
a single variable r. Three widely known representations 
for r were used for the fit. The experimental data for pi- 
ons on nitrogen and for identified hadrons on krypton 
nuclei obtained by the HERMES experiment [3, 4] were 
used for this fit. We have demonstrated that the nuclear 
attenuation ratio can be presented, with good precision, 
as a function of a single variable t instead of a function 
of two variables v and z. Moreover, i?^^ is a linear func- 
tion of r. We named t a scaling variable, because it con- 
tains all V and z dependencies of For the fit we 
have obtained as a function of r from experimen- 
tal data UrU], where it was measured as a function of v 
with the integration over z, and as a function of z with 
the integration over v. This means, that the data were 
taken at an unequal binning over v and z. In case of in- 
dependence the detailed bins over the variable v were 
taken, and for each value of v the value of z averaged 
over whole range of measured z {< z >) was taken. In 
case of z-dependence we also performed the above men- 
tioned procedure. 

Recent work of published by the HERMES experi- 
ment [6] allows to escape this difficulty, because the data 
published contains, among others, the so called two- 
dimensional data, i.e. multiplicity ratio i?5vf for charged 
pions produced in a z slices as a function of v, and in a 
V slices as a function of z. The data were obtained for 
four nuclear targets and used for a new fit of R\i as a 
function of t. 

The main aim of this work is to show, that R\[ is a 



1 In fact, R'^j also depends on ttie photon virtuality and on the 
square of ttie hadron transverse momentum in respect to the virtual 
photon direction, p^. However, from the experimental data, it is 
known that Bl^j is a much sensitive function of v and z in compari- 
son with and pf. 
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function of single variable r, rather than a function of 
two variables v and z, using the new set of given 
by the two-dimensional analysis, where the data is split 
into more regular v and z bins than in case of traditional 
presentation in form o{ v - and z - dependencies. This 
will allow to verify the results of our preceding work 1^ 
in more favourable conditions, i.e. to confirm that in 
electroproduction of hadrons in nuclear medium we in- 
deed observe scaling, where r takes over the role of the 
scaling variable. 

This paper is organized as follows. Nuclear atten- 
uation in an absorption model is presented in the next 
section. In Section 3 we discuss the choice of an appro- 
priate form for the variable r. Section 4 presents results 
of the fit. Conclusions are given in Section 5. 



II. NUCLEAR ATTENUATION IN ABSORPTION 
MODEL 

The semi-inclusive DIS of lepton on nucleus of 
atomic mass number A is 



(1) 



where liilf) and h denote the initial (final) leptons and 
the hadron observed in the final state. This process is 
usually investigated in terms of : 



B!1^{v,z) = 2daA{v-,z)IAdGo{v,z)- 



(2) 



Experimental data are usually presented at precise val- 
ues of one variable and average values of another^ . In 
this work we adopt a model according to which the ori- 
gin of the nuclear attenuation is the absorption of the 
prehadron (string, dipole) and final hadron in the nu- 
clear medium. In that case E!\^ has the following form: 



p{b,x)[Wib,x)]'^'^-^Ux, 



(3) 



where b is the impact parameter and x the longitudi- 
nal coordinate of the DIS point, p is the nuclear density 
function with a normalization condition j p{r)d^r = 1. 
W{b, x) is the probability that neither the prehadron nor 
the final hadron h are absorbed by a nucleon located 
anjrwhere in the nucleus. For W{b, x) we use the one- 
scale model proposed in Ref. 10]: 



(4) 



W{b, x)^l-aq Pq{x' - x)p(b, x')dx' 

J X 

-Oh / Ph{x' - x)p{b,x')dx' , 



where (iq and ah are the inelastic cross sections for 
prehadron-nucleon and hadron-nucleon interactions, 
respectively. Generally speaking cjq is a fimction of the 
distance x' — x, v and z However a comparison of 
simple theoretical models containing aq as a parame- 
ter P, 0], with the experimental data obtained in dif- 
ferent kinematical regions, (in particular in different do- 
mains of v), shows that the approximation aq — const 
leads to a very acceptable agreement with the data. Fur- 
ther, taking into accoimt the qualitative character of this 
work, we shall use this approximation. In the region of 
moderate energies, ah are approximately constant for all 
h. Pq{x' — x) is the probability that at distance x' — x from 
the DIS point, the particle is a prehadron and Ph [x' — x) 
is the probability that the particle is a hadron. The 
abovementioned probabilities are related via a condition 



Ph{x' -X) = 1-Pq{x' -X). 



(5) 



In analogy with the survival probability for a particle 
having lifetime r in a system where it travels a distance 
x' ~ X before decaying, Pq{x' — x) can be expressed in 
the form 



Pq{x' — x) — exp[— (x' — x)/t], 



(6) 



where r is the formation time. Substituting expressions 
for Pq{x' — x) and Ph{x' — x) in eq.(4) one obtains 



W{b,x) K 1 - ah / p{b,x )dx + T{ah - aq)p{b, x) 

J X 

wi{b,x) + T{i',z)w2{b,x). (7) 

In the framework of our suppositions about aq and ah, 
W depends on ly and z only by means of t(i^, z). 

In more detail, the formation time in the string mod- 
els can be divided in two parts (see, for instance, the 
two scale model presented in Refs. The first part 

is the constituent formation time Tc, which defines the 
time elapsed from the moment of the DIS untill the pro- 
duction of the first constituent of the final hadron. The 
second time interval begins with the production of the 
first constituent until the second one, which coincides 
with the yo-yo"* or final hadron production. Compari- 
son with the experimental data shows that in the sec- 
ond interval, the prehadron-nucleon cross section has 
values close to the hadron-nucleon cross section ah- If 
the difference between these cross-sections is neglected, 
the model is reduced to one scale model with r = Tc. 
Substituting W{b, x) in R\,j we obtain 



R 



M 



d^b 



p{b,x){'wi + T'W2)^^ ^^dx 



an + Taa + r-a^ + 



(8) 



^ In case where the u - dependence is studied, we have R\,j{i', (z)), 
where {z) are the average values of z for each u bin. For z - depen- 
dence R'^j{{u), z),where (u) are the average values of u for each z 
bin. 



^ (Tq is a function of formation time r rather than a function of vari- 
ables u and z separately. 

■* The yo-yo formation means that a colorless system with valence 
contents and quantum numbers of the final hadron is formed, but 
without its "sea" partons. 
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where i is the maximal power of t with which we are 
Hmited. Although R']^^ is a polynomial of r with maxi- 
mal power A— 1, it is expected that a^i > ai2 > > ■ ■ ■. 
The coefficients ajj depend on A, a-q, cr/j and nuclear den- 
sity. For fit we use three expressions for R'Ij as first, sec- 
ond and third order polynomials of r : 



Rm[Pi] = an +Tai2, 



Rm[P2] = 021 + ra22 + r^a23, 



(9) 



(10) 



RmIPs] = 031 + ra32 + r^ajs + T^asi- (H) 

In order to get the information on the influence of high- 
est order polynomial forms for i?^, i?^[P4] expression 
also was checked (see section 4). 



III. FORMATION TIME 



Equation (8) shows that within our approximation, 
depends on i' and z only by means of r(j/, z). This is 
the reason why we call r a scaling variable. In this sec- 
tion we shall discuss the physical meaning and possible 
expressions of the formation time t. There are differ- 
ent definitions for the formation time. We define it as a 
time scale which is necessary for the prehadron-nucleon 
cross section to reach the value of the hadron-nucleon 
cross-section. In the literature there are three qualita- 
tively different definitions for r. In the first extreme case 
it is assumed that r = (Glauber approach). In the sec- 
ond extreme case r ^ ta, where rA is the nuclear radius 
(energy loss model [9]). And at last, and in our opinion 
more realistic definition of the formation time, as a func- 
tion of u and z which can change from zero up to val- 
ues larger than ta- Experimental data seem to confirm 
that for moderate values of ly (on the order of IQGeV) 
the formation time is comparable with the nuclear size, 
i.e. the hadronization mostly takes place within the nu- 
cleus. This follows from the comparison of the expe- 
rimental data for R'Ij obtained in the region of moder- 
ate 1 4] and high |2] energies. At moderate energies R'Ij 
significantly differs from unity and is a sensitive func- 
tion of ly and z, at high energies R'^j w 1 and weakly 
depends on i' and z. For the formation time we shall use 
expressions which do not contradict the third definition 
mentioned above. The following expressions are used: 
1. Formation time for the leading hadron [ ,10.] , which 
follows from the energy-momentum conservation law 



in its rest frame tq and Lorentz factor (see, for instance, 
Ref. 17]) 



TLo 



TO- 



To- 



ruh 



(13) 



where Eh and nih are the energy and mass of the hadron 
h, respectively. In Ref. we have discussed forma- 
tion time TLor. in detail. In particular we considered the 
possibility of that tq being proportional to rn;,. In the 
present paper we have to deal with hadrons of one type 
only - charged pions - and looking at interested u and 
z dependencies of r, because all factors, which do not 
depend from these variables, can only renormalize coef- 
ficients fly (see equations (9)-(ll)). 

3. The formation time following from the Lund string 
model in Ref. UJ is^ 



'^Lund 



ln(l/z")-l + z^ 
l-z2 



zv 

K 



(14) 



One should note that all three types of formation 
time have similar behavior with v, but different behav- 
ior with z. At the values of z typical for the HERMES 
kinematics (z > 0.2) the behavior of t defined as in 
eqs.(12) and (14) with z is similar, i.e. they are decreas- 
ing with the increase of z, while t defined as in eq.(13) 
is increasing with the increase of z. 



IV. RESULTS 

The two-dimensional data from 0, i.e. the multi- 
plicity ratio R^li for charged pions produced on helium, 
neon, krypton and xenon nuclei, in a z slices as a func- 
tion of V, and in a slices as a function of z were used 
to perform the fit. The independent fit including 47 ex- 
perimental points was performed for each nucleus. As 
it is clear from eq.(2), experimental points correspond- 
ing io V - dependence R\[{v, (z)), and z - dependence 
R\j {{v),z) enter in fit on equal basis, as a values of func- 
tion R\,j{v., z) at values of variables {v, z) equal {v, (z)) 
and {{v), z), respectively. For the fit R\j has been taken 
in polynomial forms R\[[Pi,2.3\ [see eqs.(9)-(ll)], and 
formation times (lengths) as in eqs.(12)-(14). The results 
for the reduced denoted as /d.o.f. presented in Ta- 
ble HI One can see that for each choice of formation time 
and for each nucleus, the values of /d.o.f. are close for 
the polynomial approximations i?J^[Pi], R\.[[P2\, and 
R\j [P3], which means that the inclusion in consideration 
of the higher order pol3momials of r does not essentially 
improve the description of the data. 



Tlead. = (1 - z)///k, 



(12) 



where n is the string tension (string constant) with nu- 
merical value K = IGeV / fm. 

2. Formation time for the fast hadron, which is com- 
posed of characteristic formation time of the hadron h 



^ Note that this approximation is used only for the sake of conve- 
nience. For numerical calculations we use the precise expression for 
TLund following from equation ti^^^^ = Ty — zv/ k with Ty taken 
from eq.(4.21) of Ref. 
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FIG. 1: Multiplicity ratios for charged pions on different nuclei as a 
function of riead. (left panels), tlot. (central panels), TLund (right 
panels). Normalized values x — t /T{max) for all r are used. On 
panels from upper to lower results for helium, neon, krypton and 
xenon nuclei are presented, respectively. Solid, dashed and dotted 
curves are results of linear, quadratic and cubic polynomial fits. The 
numerical results for the linear fit presented on the panels. Experi- 
mental points obtained from 16]. Details see in text. 



mately proportional to v and 1 — z. In contrary, vari- 
able TLor. is proportional to v and z, and as a conse- 
quence R'ljii', (z)) and R\;[{v),z) have opposite behav- 
ior as functions of t^o,. . This means, that without any 
calculations one can state that nead. and TLund can serve 
as a scaling variables, but tlot. cannot. For the sake of 
convenience we have renormalized t to x — t /T{max), 
where T{max) are the maximum values of t for each set 
of data and each choice of the r expression. 

The range of variation of r and the numerical val- 
ues for T{max) in all scenarios for instance on krypton 
are: 2.46 < tlot. < U.Gfm, 0.37 < Tiead. < 15.5/to, 
0.36 < TLund < 8.28/m. Presentation of i?^^ as a func- 
tion of X allows us to place all data in an interval (0, 1). 
This choice does not influence the results of the fit and 
On the figure the linear polyno- 



the values of R'^l- 
mial is presented an 



*12 



with values an and 



'12 



ai2T{max) corresponding to the best fit. Solid, dashed 
and dotted curves represent the i?^^[Pi], i?J/[P2]/ and 
R'1,j[Ps] polynomial fit, respectively One can easily see 
that the difference between the curves corresponding to 
Rm[Pi], R'm[P2], and R'ljiPs] is small. The vertical posi- 
tions of the experimental points are the same in all sce- 
narios. The experimental points can be closer together 
(or not) depending on the type of the formation time 
definition. When looking at the t dependencies, the 
points corresponding to I'-dependence preserve their or- 
der in all scenarios. In case of ^-dependence, the order 
is the same in the scenario with tlot., but changes to op- 
posite in other scenarios. 

As a last remark, one should note that results of this 
analysis do not depend from the values of the parame- 
ters, in particular from the values of k and tq. 



In order to test this, we have also calculated the 
P-MlPi] polinomial form and obtained the values of 
/d.o.f. close to the ones in case of R'IjIPs]. From Ta- 
ble H] one can see that the fit gives unexpectedly good 
values for x^/d.o.f. in case of n^ad. and TLund, for tlot. 
the agreement is much worse. As it is known from ex- 
periment (3- [3-0]/ R\[(yi (z)) increases with increasing 
of u, and R\f{{v), z) decreases with increasing of z for 
all nuclei. Our assumption is that these functions in- 
deed present different representations of the same func- 
tion, which depends on variable t only. In turn t is a 
function of v and z. Now let us discuss the figure. Us- 
ing eqs.(12)-(14), we present R^^iv, (z)) and R'Ij{(i'),z) 
as fimctions of r. Experimental points and results of 
the fit are presented in Fig. [T] Solid points correspond 
to the R\i{v, (z)) obtained from the experimental data 
for i/-dependence, open points to the R\j{{v), z) from z- 
dependence. From the figure one can easily note that ex- 
perimental points for R\[ [u,{z)) and R\.j {{v),z)as func- 
tions of r have the same behavior and approximately 
coincide when rjead. and TLund serve as the variables. 
The reason for this is that these variables are approxi- 



V. CONCLUSIONS. 

The two-dimensional nuclear attenuation data for 
charged pions on helium, neon, krypton and xenon nu- 
clei obtained recently by the HERMES experiment 0] 
were used to perform the fit. So far it has been supposed 
that experimentally measured function R\i{v, (z)) de- 
pends from variable v only, and R'Ij{{i^), z) from z only 
In our preceding work we have assumed that these 
functions indeed present different representations of the 
same function, which depends on only one variable t. 
In turn t is a function of v and z. In this work based 
on two-dimensional distributions (more suitable for this 
investigation), the results of our preceding work |5] are 
confirmed. We demonstrate (see Table 1 and figure) that 
R\[{v, (z)) and z) as a functions of r have the 

same behavior and approximately coincide, when r is 
used in form of Ti^ad. and TLund- In contrary, they have 
opposite behavior as a functions of TLor. ■ This indicates 
that T/ead. and TLund Can serve as scaling variables, but 
TLor. cannot. We also show that R\.j{t), with a good pre- 
cision, can be parametrised in a form of a linear polyno- 
mial an + rai2, where the fitting parameters an and ai2 
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TABLE I: The jd.o.j. values obtained from polynomial fit. Pi, 2, 3 denote the expressions -Rm[A,2,3] used as fitting functions. 
The necessary details concerning the data sets used for the fit are given in the text. 



'^lead. '^Lor. '^Lund 



A 


exp 


Had. 


Pi 


P2 


Pa 


Pi 


P2 


^3 


Pi 


P2 


^3 


"He 


47 


< TT > 


0.43 


0.42 


0.39 


0.37 


0.37 


0.38 


0.42 


0.39 


0.36 


20iVe 


47 


< TT > 


0.62 


0.55 


0.50 


0.98 


0.98 


0.96 


0.45 


0.38 


0.35 




47 


< TT > 


1.27 


0.83 


0.76 


7.18 


7.31 


7.47 


0.73 


0.57 


0.50 




47 


< TT > 


0.83 


0.55 


0.50 


7.68 


7.85 


7.98 


0.54 


0.45 


0.42 



do not depend on ly and z. We conclude that experi- 
mentally measured function is a function of single 
variable t, and that t scaling follows naturally from the 
absorption model, but we do not confirm that scaling 



is a property of an absorption model only. There could 
exist other mechanisms which can also lead to the scal- 
ing behavior - this is an open question, which requires 
further investigation. 
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